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Abstract. A large area ring-imaging Cherenkov detector has been designed to provide clean hadron iden¬ 
tification capability in the momentum range from 3 GeV/c up to 8 GeV/c for the CLAS12 experiment at 
the upgraded 12 GeV continuous electron beam accelerator facility of Jefferson Laboratory. The adopted 
solution foresees a novel hybrid optics design based on aerogel radiator, composite mirrors and highly 
packed and highly segmented photon detectors. Cherenkov light will either be imaged directly (forward 
tracks) or after two mirror reflections (large angle tracks). We report here the results of the tests of a large 
scale prototype of the RICH detector performed with the hadron beam of the CERN T9 experimental 
hall for the direct detection configuration. The tests demonstrated that the proposed design provides the 
required pion-to-kaon rejection factor of 1:500 in the whole momentum range. 


PACS. 29.40.-n Photomultipliers in nuclear physics - 
Cherenkov detectors 

1 Introduction 

The Jefferson Laboratory (JLab, USA) main facility is 
currently completing a major upgrade program that al¬ 
ready doubled the energy of the electron beam (from 6 
GeV to 12 GeV). This program will lead in the near fu¬ 
ture to an increased luminosity, the enhancement of the 
detector systems in the three existing Halls A, B and G 
and the construction of a new experimental Hall D. The 
Hall B will house the new CLAS12 large acceptance multi- 
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purpose spectrometer, that will operate with a highly po¬ 
larized electron beam with energy up to 11 GeV and a 
luminosity as high as 10^^cm“^s“^ impinging on unpolar¬ 
ized as well as polarized hydrogen and deuterium targets, 
thus providing unique conditions for the study of electron- 
nucleon scattering [T]. 

The CLAS12 experimental program [2] covers many 
topics in hadron physics. Particular attention is devoted 
to the 3D imaging of the nucleon through the study of 
generalized and transverse momentum dependent parton 
distributions (GPDs and TMDs) in the poorly explored 
valence region (high Bjorken x). A number of approved 
experiments [3] require an efficient hadron identification 
in the 3-8 GeV/c momentum range in order to allow the 
flavor separation. As the flux of pious is one order of mag¬ 
nitude larger than kaons, a pion rejection factor of about 
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1 : 500 is required to limit the pion contamination in the 
kaon sample to a few percent level. 

The standard particle identification (PID) is performed 
in CLAS12 by using a time-of-flight (TOP) system and 
two Cherenkov gas detectors with high (HTCC) and low 
(LTCC) thresholds. The TOF is able to provide full hadron 
PID up to about 3 GeV/c, while the two Cherenkov de¬ 
tectors reach the required pion rejection factor only at the 
edge of the available phase space (hadron momenta above 
7 GeV/c) and are not able to distinguish kaons from pro¬ 
tons. In order to achieve the neccessary PID requirements 
in the full momentum range, a modular Ring Imaging Che¬ 
renkov (RICH) detector has been proposed to replace at 
least two of the six azimuthal sectors of the LTCC. Each 
RICH sector will have a projective geometry, with a depth 
of about 1.2 m and almost 5 m^ entrance window area. 

Based on simulation studies EE], the best configu¬ 
ration for the detector is a non-conventional proximity- 
focusing design, with a wall of aerogel radiator tiles, an 
array of visible light photon detectors and a mirror sys¬ 
tem. The latter is essential to reduce to about 1 m^ per 
RICH sector the area covered by the photon detectors, 
thus minimizing costs and the material-budget impact on 
the downstream CLAS12 detectors (TOF and Calorime¬ 
ters). 

The concept of this RICH is illustrated in Fig. [H For 
forward scattered particles {0 < 12^) with momenta in the 
range between 3 and 8 GeV/c, the Cherenkov light will be 
directly detected by the photon detector array (Fig.[T]top). 
For particles with larger incident angles (12^ < 0 < 35*^), 
with momenta between 3 and 6 GeV/c, the Cherenkov 
light will be double-reflected by a spherical and a planar 
mirror and focused onto the photon detector array after 
two passages through the lower section of the aerogel wall 
(Fig. [T] bottom). This double pass imposes the use of aero¬ 
gel with reduced thickness (2 cm) and with relatively high 
refractive index (n=1.05) in the lower section of the detec¬ 
tor, to minimize the photon yield loss. The photon yield 
losses will be compensated by the use of a thicker (6 cm) 
aerogel in the upper section. 

A variety of laboratory and beam tests have been per¬ 
formed in order to validate this innovative configuration. 
Laboratory tests have been done to study the single RICH 
component performances, in particular the optical proper¬ 
ties of aerogel produced by various companies worldwide 
and the performaces of different high-granularity photon 
detectors in the Single Photo-Electron (SPE) regime. In 
addition, a campaign of beam tests has been pursued using 
prototypes of increasing scale and complexity, performed 
at the Frascati Beam Test Facility [7] and at the CERN T9 
beam line [8]. In the following, we report the results of the 
final tests performed at CERN for the direct light configu¬ 
ration, including several tests performed on the main com¬ 
ponents, whose performance was challenged in the limit of 
high momenta. The prototype used in these tests was de¬ 
signed in such a way that it reproduces the direct light 
configuration of the CLAS12 RICH. 

For the reflected light configuration, the prototype was 
designed following the same concept of the CLAS12 RICH 



aerogel 



aerogel 

Fig. 1. The conceptual design of the RICH detector for for¬ 
ward (top plot) and large angle (bottom plot) particles. 


(i.e. the use of spherical and planar mirrors and the double 
pass through the aerogel), although geometry limitations 
imposed a reduced light path length, not optimized for real 
hadron separation. Preliminary results for the reflected 
light configuration have been already published [6], the 
final results will be reported in a future paper. 


2 The experimental setup 

2.1 The CERN T9 beam line 

The setup for the tests was installed at the T9 experimen¬ 
tal area located in the East Area of the PS/SPS complex 
at CERN [8]. The primary proton beam of the PS is sent 
to a target to produce a secondary beam. Different targets, 
magnets and collimators allowed to select the composition 
(electrons, muons or hadrons), the charge and the momen¬ 
tum [9] of the secondary beam. 

The tests were performed using the settings for nega¬ 
tive charge hadrons with the highest momentum possible, 
allowing to cover the region between 6 and 8 GeV/c, the 
most critical one for the CLAS12 RICH. The relative pop¬ 
ulation of : K~ : p was approximately 160 : 5 : 1, ba¬ 
sically independent from the momentum in the range be- 
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tween 6 and 8 GeV/c. A number of dipole and quadrupole 
magnets with predetermined settings allowed to focus the 
beam at different distances from the beam pipe exit win¬ 
dows. The choice was to have the smallest beam spot as 
close as possible to the RICH photodetector plane, at 10 
meters from the exit window of the beam line. 

The time structure of the beam was determined by the 
number of spills (from 1 to 3) extracted from the primary 
PS proton beam. Each spill had a duration of 400 ms and 
the period of the PS operation was 40 s, corresponding to 
an overall duty cycle of the order of few percent. 

The T9 equipment included a CO 2 gas Cherenkov coun¬ 
ter used to separate pious from heavier hadrons by ad¬ 
justing the gas pressure so that only pions were above 
threshold. 


2.2 The RICH prototype 

The RICH prototype has been designed and built at the 
INFN Laboratory of Frascati. It comprises a large (ap¬ 
proximately 1.6 X 1.8 X 1.6 m^) light-tight box and in¬ 
ternal modular supports holding the various components, 
that may be adjusted depending on the measurement con¬ 
ditions. The proximity gap length was about 1 m, close to 
the CLAS12 conditions. A schematic drawing of the setup 
and a picture of the RICH box during the installation are 
shown in Fig. [2j 


2.2.1 The radiator 

An excellent radiator for RICH hadron identification in 
the few GeV momentum range is silica aerogel, an amor¬ 
phous solid network of Si 02 nanocrystals with a very low 
macroscopic density and a refractive index in between 
gases and liquids. It has been used as radiator material 
for RICH detectors in several particle physics experiments 
[iniiiiiiniiia and is planned for future use El- 

A systematic characterization has been carried out on 
a variety of aerogel samples from different manufactur¬ 
ers. The aerogel from the Budker and Boreskov Cataly¬ 
sis Institutes of Novosibirsk p!^[T6] . that combines high 
transparency with a considerable geometrical (area and 
thickness) flexibility and mass production capabilities, was 
chosen. Noteworthy, during the prototyping, the produc¬ 
tion technique and the resulting quality have significantly 
improved in time, reaching a clarity parameter as low as 
0.0050 /im^ cm“^ for a n = 1.05 refractive index at the 
reference wavelength Aq = 400 nm (blue light). 

Precise measurements of the aerogel transmittance have 
been performed, using a Lambda 650 S PerkinElmer spec¬ 
trophotometer, on a broad range of light wavelengths be¬ 
tween 200 and 900 nm [T7]. Typical transparency values 
in the blue light region are close to 70% for the most re¬ 
cent production of tiles with refractive index of 1.05 and 
20 mm thickness, while lower values, close to 60%, have 
been measured for the older production of aerogel with 
refractive index of 1.04 and 1.06. 


GEM1 


"~W~ 

Trigger 

scintillators 





Beamline 
exit window 



Fig. 2. Top plot: a schematic drawing of the setup in the 
T9 experimental hall. Bottom plot: picture of the RICH box 
installed at the CERN-T9 beam line; the beam enters from the 
right; the photon detector array is on the left of the picture, 
while the support for the aerogel tiles is on the right. 


The Cherenkov light emission spectrum per unit radi¬ 
ator thickness x and wavelength A can be written as 


d^Nc 

dxdX 


2TTa 




( 1 ) 


where rjc is the Cherenkov emission angle. The spectrum 
has a broad distribution and it is rapidly increasing when 
the wavelength changes from the visible light toward the 
UV region. Thus, it is important to keep under control the 
chromatic dispersion, i.e. the dependence of the refractive 
index on the light wavelength, as it constitutes one of the 
largest contributions to the Cherenkov angle resolution. 
One of the standard techniques for its measurement is 
the so-called prism method in conjunction with the Snell- 
Descartes formula, but experimental measurements exist 
in literature only for refractive index up to n = 1.03 [18]. 

In our prototype, we used aerogel tiles with an area of 
approximately 58 x 58 mm^, 20 mm thickness and refrac¬ 
tive indices n = 1.04, 1.05, 1.06. The Novosibirsk aerogel 
is hydrophilic, i.e. tends to absorb water molecules from 
the air, resulting in a degraded transparency. Thus, it was 
important to monitor the transparency of each individual 
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aerogel tile during the data taking and eventually correct 
the data for the degraded transparency. For this reason, a 
strict protocol to store and handle the aerogel was estab¬ 
lished. In addition, before and after the installation, the 
relative transparency of each tile was monitored using a 
simple setup made by a blue light laser (A = 405 nm), 
a photodiode and a multimeter. The transmittance was 
then evaluated as the ratio between the multimeter read¬ 
out with and without the aerogel. This method has the 
advantage of being very fast (a few minutes for each set of 
measurements) but provides only relative values of trans¬ 
parency. The results showed only marginal degradation in 
the transparency between the measurements performed 
before and after the usage of the tiles, thus no corrections 
have been applied to the data. 

2.2.2 The photodetectors 

Dedicated simulation studies [?] have shown that the re¬ 
quired Cherenkov angle resolution can be achieved if the 
pixel size of the photon detectors is less than 1 cm. 

The Hamamatsu H8500 multianode photomultiplier 
tubes (MAPMTs) [19] have been selected as a good candi¬ 
date for the first module of the CLAS12 RICH, being an 
effective compromise between detector performance and 
cost. It comprises an 8 x 8 array of pixels, each with di¬ 
mensions 5.8 X 5.8 mm^, into an active area of 49 x 49 mm^ 
with a very high packing fraction of Cpack ~ 89% and a 
collection efficiency Ccoii ~ 80% [20]. The device offers a 
spectral response matching the spectrum of light trans¬ 
mitted by the aerogel, with a quantum efficiency peaking 
at 400 nm, and a fast response (less than 1 ns rise time) 
useful to suppress the background. 

Although the H8500 MAPMT is not advertised as the 
optimal device for single photon detection purposes, the 
characterizations performed in laboratory on several units 
in this regime [21] have demonstrated that this device 
can achieve performances adequate for the CLAS12 RICH 
requirements. The uniformity of the H8500 response has 
been extensively studied with a pico-second pulsed laser. 
The typical gain variations in the pixel response of each 
MAPMT, of the order of 1:2, can be easily compensated 
by the readout electronics. Sub-millimeter precision scans 
have been used to study the PMT response across the 
MAPMT surface and to evaluate the true active areas 
of the pixels [22]. The fraction of SPE signal below the 
pedestal threshold has been reduced to less than 15% op¬ 
erating the high voltage at 1040V or above, while less than 
5% crosstalk has been measured. 

The RICH prototype was instrumented with 28 MA¬ 
PMTs, mounted on a circular support, that can be radi¬ 
ally moved so to intercept the Cherenkov rings produced 
with different opening angles depending on the aerogel 
refractive index. Fourteen Hamamatsu H8500C, with nor¬ 
mal borosilicate glass window (wavelength cut-off at about 
250 nm), and fourteen Hamamatsu H8500C-03, with glass 
window of enhanced transparency for UV photons (wave¬ 
length cut-off at about 185 nm) were alternated along 
the ring. The ring coverage varied between about 80% 


for n = 1.04 aerogel refractive index to about 60% for 
n = 1.06. 

In order to increase the separation of the SPE region 
from the pedestal, the high voltage of the MAPMTs was 
set at 1075V, close to the recommended upper limit of 
llOOV. 


2.3 The tracking system 

Particles’ trajectories were reconstructed using two 10 x 10 
cm^ triple-GEM chambers [23]. The GEM foils and the 
double layer readout strip plane were produced at CERN 
M- Each chamber measures the particle hits on a (x, y) 
plane made of 256 x 256 strips, with a nominal spatial 
resolution of about 100 ym. 

The two chambers were placed about 4.2 m apart from 
each other, with the upstream one placed 3.5 m before the 
MAPMT plane and the downstream one 0.7 m behind the 
MAPMT plane. 


2.4 The trigger 

The trigger of the Data Aquisition (DAQ) system was 
produced by two small plastic scintillators, partially over¬ 
lapping on a small spot of about 2x3 cm^, placed 1 m 
downstream of the RIGH MAPMT plane and just after 
the second GEM chamber. The scintillators were readout 
by Hamamatsu R1450 single anode PMTs. The output 
signals were amplified, discriminated and then put in co¬ 
incidence with the signal coming from the beam by using 
standard NIM electronics. This coincidence was then trig¬ 
gering the readout electronics of the various systems. 


2.5 Readout electronics 

The prototype was equipped with 28 MAPMTs with a 
total of 1792 channels, thus requiring highly integrated 
and compact readout electronics. We used custom-made 
electronics based on the 64 channel MAROG3 [25] chip, an 
ASIG manifactured by the Omega Group of LAL (Paris, 
France) for the ATLAS experiment. 

The chip embeds an individual channel preamplifier 
stage, with variable gain, followed by a highly configurable 
shaping section, and provides binary and analog charge 
measurements on independent lines. The binary output 
is produced with an adjustable threshold discriminator 
array, 64-channels wide, while the analog is obtained with 
a sample and hold circuit and it is serially multiplexed out 
and digitized by an external ADC. During this test, only 
the serial analog output has been used. 

Each MAROC3 chip is mounted on a front end card 
connected to the MAPMT through high density coaxial 
fiat cables and equipped with an FPGA in charge of the 
configuration, the control and the readout of the chip. 
The 28 front end cards are installed on 4 backplanes that 
provide electrical connectivity for power supply and data 
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Fig. 3. Top plot: the various components of the readout elec¬ 
tronics. Bottom plot: the assembled setup installed in the 
RICH prototype box. 


transfer. Two control boards act as event builder and pro¬ 
vide USB2.0 interface to the aqnisition node. The total 
number of available channels is 4096 in a very compact 
system [26]. The various components and the assembled 
set-up are shown in Fig. [3l 

The GEM electronics [27] is based on dedicated VME 
modules and on the analog APV25 chip developed by 
Imperial College for the CMS silicon detectors [28] and 
first used in a GEM detector by COMPASS. The APV25 
is designed to be radiation tolerant (at least up to 10 
Mrad), has a high channel density (128 channels/chip) 
and a rather fast readout (up to 40 MHz). Each APV25 
channel consists of a preamplifier followed by a shaper and 



Fig. 4. Typical ADC spectrum of the threshold Cherenkov 
counter. The dashed red line indicates the threshold to separate 
pions from heavier hadrons (5 a above the pedestal mean). 


a 192 cell analog pipeline, where the shaped signal is con¬ 
tinuously sampled at 20 or 40 MHz. The samples awaiting 
readout are flagged by external triggers, with adjustable 
latency. Multiple samples (up to 6) per channel and event 
will be converted and acquired by the GEM electronics. 
The acquired information permits a simple signal analysis 
which provides time (with resolution of about 5 ns) and 
charge information. 

The DAQ hardware electronics was completed by a 
CAEN V785N ADC for the threshold Cherenkov counter 
readout. 

The DAQ software managed the acquisition of the dif¬ 
ferent systems through a CAEN V2718 VME bridge, con¬ 
nected to the computer via optical link. Events for each 
electronic system are flagged by the DAQ with a common 
event number which allows the full and unambiguous re¬ 
construction of the event. 


3 Data reconstruction 

3.1 The gas threshold Cherenkov counter 

A typical spectrum of the gas threshold Cherenkov counter 
is shown in Eig. H] where the broad distribution is due to 
the pions while the narrow pedestal peak is due to hadrons 
below the Cherenkov emission threshold. The cut at 150 
channels (5cr above the pedestal peak) is set to separate 
pions from heavier hadrons. 

3.2 The photodetectors 

The electronic system based on the MAROC3 chip al¬ 
lowed to individually record ADC information of all the 
28 X 64 = 1792 pixels of the MAPMTs. In order to min¬ 
imize the SPE loss, a uniform gain ^ = 4 has been set 
to all the MAROC3 channels. No finer channel-to-channel 
equalization has been performed. 

A typical ADC spectrum registered from one of the 
channels is shown in Eig. [S] the pedestal and the SPE 
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Fig. 5. Typical ADC spectrum (in logarithmic scale) of one 
of the MAPMT pixels registered by the MAROC3 electronics. 
The histogram has been fitted with the sum of two gaussian 
curves for the pedestal and the first photoelectron peak, shown 
by the curves in the plot. The fitted parameters are reported 
in the legend and described in the text. 


regions are clearly visible and well separated. The figure 
shows also a fit with the sum of two gaussians, one for 
the pedestal and one for the SPE signal, weighted by the 
Poisson probability to have zero or one photoelectron, re¬ 
spectively. In the legend, the values of the gaussian mean 
and sigma of the pedestal {Qq and Sq) and of the SPE 
peak (Qi and Si) as well as the average number of photo¬ 
electrons /i and the normalization constant are reported. 
A cut at 5cr above the pedestal peak selects the Cherenkov 
photon hits. The SPE efficiency has been estimated by ex¬ 
trapolating the SPE gaussian of the fit in Eig. 0 below the 
pedestal cut. We obtained espE ~ 0.87. 

The ADC spectra of all the pixels have been carefully 
analysed in order to identify and eventually remove noisy 
or dead channels. At the end of this data quality analysis, 
about 20 bad channels have been found, corrresponding 
to about 1% of the total. 

Eor each run, an event display allowed to monitor on¬ 
line the quality of the data. Examples of the measured hit 
distributions for runs with different aerogel refractive in¬ 
dex n are shown in Eig. [6l Rings of increasing radius from 
n = 1.04 to n = 1.06 can be seen. 
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3.3 The GEMs 

The GEM tracking system consists of a pair of chambers, 
each one with readout on the (x, y) strip planes, located 
upstream and downstream of the RICH. Eor each strip, 
three samples of the charge signal, separated by 25 ns, are 
recorded when a trigger is issued. The track reconstruction 
proceeds through the following two steps. 

The first step is the identification, in a given GEM 
plane, of a peak at the same location (i.e. same strip) in 
all the three time samples. A typical example of one event 


Fig. 6. Online monitoring of the Cherenkov photon hit pat¬ 
terns measured with aerogel of different refractive index n. 
Data with 8 GeV/c beam momentum. 
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Fig. 7. Typical GEM event readout on one plane: ADC signal 
as a function of the strip ID, after pedestal subtraction, at the 
three sampling times Tl, T2=Tl+25 ns and T3=Tl+50 ns. 



Fig. 8. Top plots: beam profiles measured on GEMO (upstream 
to the RICH) in the horizontal (left plot) and vertical (right 
plot) directions for several runs. Bottom plots: the same for 
the GEMl (downstream to the RICH). Data for 8 GeV/c mo¬ 
mentum. 


for one GEM plane is shown in Fig. [71 In this case, a good 
signal is easily identified on the strip 169 in all the three 
samples. All the peaks exceeding 100 ADC units above 
the pedestal are ordered according to their amplitude and 
kept for further analysis. The second step is the matching 
of the X and y peaks in one chamber to obtain a good 
GEM hit and finally to match a pair of good hits on the 
two GEM chambers to reconstruct a good hadron track. 

The beam profiles measured in the two GEM cham¬ 
bers are shown in Fig. [8] for several different runs over the 
data taking period. No major variations in the beam pro¬ 
file have been found during the test period, ensuring the 
stability of the beam conditions and of the tracking sys¬ 
tem. The holes in the distributions are due to noisy strips 
that have been removed in the off-line data quality analy¬ 
sis. The profile is wider on GEMO (upstream GEM) than 
on GEMl (downstream GEM) because of the focusing of 
the beam. The horizontal width is about twice the vertical 
one due to a different beam collimator in x and y. 


Due to the beam structure, on average no more than 
one good charged track per trigger was expected. The frac¬ 
tion of events with more than one reconstructed track is 
about 3%. This is due to either electronic noise on the 
GEM readout induced by the MAROG3 electronics or fake 
peaks produced by the reconstruction algorithm in corre- 
spondance of the removed noisy strips. These additional 
hits have usually much smaller amplitude than the physi¬ 
cal ones. For this reason, only the track resulting from the 
hits with the biggest amplitude in each event is considered 
as a good track hit candidate. 

The hit reconstruction efficiency is above 80% on GEMO 
and around 60% on GEMl, the latter lower value due to a 
leak of gas that was not possible to completely fix. Overall, 
the track reconstruction efficiency was about 50%, lower 
than the typical value (about 80%, with a 90% efficiency 
on each chamber) expected for this type of detector. 


4 Event reconstruction 

4.1 Alignment of the tracking system 

Rough alignment of the GEM chambers was done during 
the installation of the setup. A finer alignment has been 
obtained off-line by fitting the GEM tracks to the ring 
reconstruction. 

A set of 4 alignment constants (Xq,Iq, k^i) has 
been introduced to take into account possible displace¬ 
ments of the GEM on the transverse plane with respect 
to the beam line. These constants have been determined 
by comparing the projection of the GEM tracks onto the 
MAPMT plane with the reconstruction of the Gherenkov 
ring performed using the MAPMT information only. 

To do this, the Gherenkov rings are reconstructed by 
minimizing the quantity 

^pe 

S{R,Xc,Yc) = J2[{xi - Xcf + {vi - Ycf - R^f ( 2 ) 

where (xi^yi) are the coordinates of the ith photon hit, 
{Xc^Yc) are the coordinate of the ring center and R its 
radius. The event-by-event minimization of eq. m with 
MINUIT [29] can be performed either to determine the 
three free parameters Xc, Yc and R {Spar fit) or to fit 
the radius only {Ipar fit) fixing Xc and Yc by using the 
GEM information. 

The four GEM alignment constants have been deter¬ 
mined by minimizing the quantity: 

Q(Vo,ro, Vi,n) = 4 + (3) 

where is the RMS of the distribution of the Gherenkov 
ring radius in the Ipar fit and AX and AY are the mean 
values of the distributions of the difference between the 
GEM track (projected on the MAPMT plane) and the 
ring center coordinates obtained in the Spar fit. 

In Fig. (9] we show an example of the Gherenkov ring 
radius distribution from the Spar fit (top, aR = 2.1 mm), 
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Fig. 9. Cherenkov ring radius from the Spar fit (top), from the 
Ipar fit before (center) and after (bottom) the GEM alignment 
procedure. Data for 8 GeV/c pions on aerogel with n = 1.05 
and t = 2 cm. 


from the Ipar fit before (center, cfr = 3.2 mm) and after 
(bottom, (Jr = 1.8 mm) the offline alignment. An improve¬ 
ment of more than 10% in the resolution has been obtained 
from the initial Spar to the final Ipar fit at the end of the 
alignment procedure. The stability of the procedure has 
been checked by calculating the alignment constants from 
several runs with different beam energies across the data 
taking period. No run dependence of these constants has 
been found. 

After the alignment, the beam spot size projected onto 
the MAPMT plane is shown in Fig. [TOl The different size 
between the horizontal and the vertical dimensions (15 
mm vertically and 25 mm horizontally) is due to the dif¬ 
ferent collimation. The small gap at x ~ — 8 mm comes 
from the noisy strips removed from the downstream GEM. 


4.2 MAPMT position correction 

For each aerogel refractive index, the position of the MA- 
PMTs on the support was moved according to the ex- 



Fig. 10. Beam spot measured with the two GEMs and pro¬ 
jected onto the RICH photon detector plane. Data for 8 GeV/c 
beam. 


pected Gherenkov ring radius and was periodically mea¬ 
sured during the data taking. Finer alignment of the MA- 
PMTs position was performed off-line by looking at the 
residuals of the photon hits, i.e. the distance of each hit 
from the fitted Gherenkov ring. In the top plot of Fig. 
m we show an example of the residual distributions for 
the 28 MAPMTs. For each MAPMT, the shift from zero 
of the mean residual represents the correction to be ap¬ 
plied to the measured position. Typically, this correction 
amounts to a few mm and never exceeds 6 mm, the size 
of the MAPMT pixel. The distributions of the residuals 
after the alignment procedure are shown in the bottom 
plot of Fig. [TTl The different widths of the residual distri¬ 
butions between the even (H8500G-03) and odd (H8500G) 
MAPMTS is due to the UV photons, as will be discussed 
in Section [5Ml 

The alignment corrections have been computed for each 
group of runs taken in the same experimental conditions. 
On average, the MAPMT position correction provides an 
improvement of the Gherenkov angle resolution of a few 
percent. 


4.3 Cherenkov event reconstruction 

The reconstruction of the Gherenkov angle has been per¬ 
formed following the method developed for the proximity 
RIGH detectors of the ALIGE HMPID [30] and the JLab 
Hall A [31] RIGHs. The method, illustrated in Fig. [T2j 
exploits the information of the tracking system and of the 
MAPMT array and provides, for a given photon hit, the 
Gherenkov angle, correctly taking into account the small 
slope of the hadron track, of the order of few mrad, which 
produces a non circular hit distribution on the MAPMT 
plane. 

The reference frame has the 2 ; axis along the beam 
line, the y axis is vertical and pointing upwards and then 
X = y X z. The thickness of the radiator is trad^ while 
tgap = 994 mm and Uair = 1.000273 are the length and 
refractive index of the air gap between the radiator and 
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Fig. 11. Hit residual distributions for the 28 MAPMTs before 
(top plot) and after (bottom plot) the position correction. 


the photodetectors. The track enters the radiator at the 
point Po and its projection onto the MAPMT plane is Pp. 
Its polar and azimuthal angles are Op and 0^, respectively. 

The Cherenkov photon is emitted at the point Pc, at a 
depth L inside the radiator, and is detected at the point P 
after the refraction when exiting from the aerogel. In prin¬ 
ciple, a recursive calculation would allow a rather precise 
determination of L. However, given the small thickness 
of the radiators with respect to the total gap length, we 
assumed L = trad I 

The point Ppcp is the projection of the emission point 
onto the MAPMT plane. The polar and azimuthal angles 
of the photon hit are 0 and 0, while the Cherenkov angle, 
i.e. the angle between the photon direction and the track, 
is 77c. It is important to note that the refraction in the 
passage from the aerogel to the air changes its polar angle 
but not the azimuthal angle. As a result, the photon is 
moving on a plane and, by using the Snell equation, if the 
angles are known one can calculate the ring radius R as 



reference frame 



1 


Fig. 12. Illustration of the Cherenkov angle reconstruction. 


The problem is solved by minimizing the difference 

g{e) = {R-f{e)f (5) 

Once 0 is obtained, the Cherenkov angle of the photon hit 
is given by 

cos(77c) = sin(^p) sin(^) cos(0 — (j)p) + cos(^p) cos(^) (6) 

Finally, the Cherenkov angle of the event, where m photon 
hits have been detected, is given by the average of the m 
Cherenkov angles. 

Knowing the Cherenkov angle, one can derive also the 
measured refractive index of the radiator 

/3cos(7?c) 

where the velocity of the particle [3 is computed from the 
beam momentum and using the gas Cherenkov counter 
information to discriminate between pions and heavier 
hadrons. 

An iterative algorithm of the Cherenkov ring recon¬ 
struction has been implemented in order to study back¬ 
ground hits. A hit is considered as background if its dis¬ 
tance from the fitted ring is bigger than a given threshold 
value. In this case, the hit is removed and a new fit is 
performed taking into account only the remaining good 
Cherenkov hits. The threshold value has been optimized 
to the data: a too loose value will include too many back¬ 
ground hits, thus worsening the resolution, while a too 
high value will remove good photons, again worsening the 
resolution. An optimal value of 12 mm (i.e. twice the pixel 
size) has been found. 


R = {Uad - L) tan(6>) + tgap - = f{0) 

sin^(6>) 

(4) 

When the emission angle is not known but the radius R 
is measured, this equation must be inverted to obtain 0. 


5 Results 

5.1 Cherenkov angle resolution studies with pions 

We discuss now the performances of the prototype in the 
nominal condition, i.e. pions with momentum of 8 GeV/c 

































10 


S. Anefalos Pereira et aL: Test of the CLAS12 RICH prototype 



Fig. 13. Total number of detected hits (dashed black line), 
the number of p.e. (full blue line) and the number of back¬ 
ground hits (dotted red line) per event measured in the RICH 
prototype with 8 GeV/c momentum pions and aerogel with 
n = 1.05 and trad = 20 mm. 


and an aerogel radiator with n = 1.05 and trad = 20 mm 
thickness. These are expected to be the most challeng¬ 
ing conditions in the CLAS12 RICH. Pions are selected 
through the gas threshold Cherenkov signal, as shown in 
Fig. m We also require a minimum of four hits on the 
MAPMTs. 

In Fig. [131 we show the total number of detected hits 
(dashed black line), the number Npe of photo-electrons 
(p.e.) for the reconstructed rings (full blue line) and the 
number Ni^kg of background hits (dotted red line). We 
have Npe = 12.02 ± 0.02 p.e. per pion Cherenkov ring and 
Nbkg ~ 1- The ring coverage of our RICH prototype was 
in this case about 70%, thus for a full coverage detection 
system we may estimate Npe ~ 17. The Cherenkov angle 
distribution is shown in Fig. [HI The gaussian fit of the dis¬ 
tribution yields rje = 309.14 ± 0.01 mrad and a resolution 
ap = 1.39 ± 0.01 mrad. 

As expected, no variations have been found in r]e as a 
function of the number of photoelectrons: the slope of the 
Npe dependence is smaller than 10“^ mrad/A/pe- On the 
other hand, the gaussian width of the Cherenkov angle 
distribution is expected to be proportional to the inverse 
of the squared root of Npe'. 


_ ^Ipe 

\^^pe 


( 8 ) 


where aipe is the single photon detection resolution. The 
Cherenkov angle resolution as a function of the number 
of p.e. is shown in Fig. [151 We measured up to 25 p.e. 
per ring and the fit with eq. ([8]) gives the single photon 
resolution aipe = 4.58 ± 0.02 mrad. 

As a test of the reconstruction algorithm, the Npe has 
been repeated by adding to eq. m a constant term. The 
results showed a negligible constant term (below 0.1 mrad) 
and basically no variation in the aipe. Thus, we can con¬ 
clude that the possible experimental artifacts, as for ex¬ 
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Fig. 14. Cherenkov angle distribution for 8 GeV/c momentum 
pions and aerogel with n = 1.05 and trad = 20 mm. The dashed 
line is a gaussian fit of the distribution, whose fitted parameters 
are shown in the table. 
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Fig. 15. Angular resolution as a function of the number of p.e. 
for 8 GeV/c momentum pions and aerogel with n — 1.05 and 
trad = 20 mm. The line is a fit with eq. with aipe = Po- 


ample residual misalignment of the setup, are well under 
control. 

The pion reconstruction efficiency may be estimated 
by computing the ratio 


atRICH 

TV 


where N^^ is the number of events with a gas Cherenkov 
counter signal above the threshold and represents 

the number of these events that are also reconstructed as 
pions in the RICH prototype. An event in the RICH proto¬ 
type is taken as a pion if the reconstructed Cherenkov an¬ 
gle falls in the range We obtain = 98.8 ±0.1%, 

with the small fraction of event loss mainly due to the 
events with less than 4 photoelectrons. 
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Fig. 16. Pion Cherenkov angle distributions (dashed black 
histograms) compared with those from events with gas Che¬ 
renkov signal below threshold (full red histograms, arbitrarily 
rescaled), for P = 6,7,8 GeV/c beam (from top to bottom). 
Data for aerogel with n = 1.05 and trad = 20 mm. 


5.2 The TT / K comparison 


Events below the gas Cherenkov counter threshold include 
both kaons and antiprotons. In Fig. [161 we compare the 
Cherenkov angle distributions measured by the RICH pro¬ 
totype for these events (full red histograms) with the dis¬ 
tributions for pions (dashed black histograms) at the three 
beam momenta P = 6,7,8 GeV/c, from top to bottom. 
The full red histograms are arbitrarily rescaled to make 
them visible. We see the prominent kaon peak separated 
from the pion one and, at the highest energy, some tail 
of the p peak, which disappears as the energy decreases 
because it goes out of the MAPMTs radial coverage. 



Hadron Momontum (GeV/c) 


Fig. 17. Top plot: mean value of the Cherenkov angle mea¬ 
sured for pions (full red circles) and kaons (empty blue squares) 
as a function of the hadron momentum, compared with the 
corresponding expected dependencies (dashed and full lines, 
respectively). Bottom plot: the same for the angular resolu¬ 
tion. Data for aerogel with n = 1.05 and trad — 20 mm. 


Gaussian fits of these distributions provide the Che¬ 
renkov angle and resolution for the different hadrons. In 
Fig- 113 we show the gaussian mean (top plot) and sigma 
(bottom plot) of the Cherenkov angle distributions as a 
function of the hadron momentum, for pions (full red cir¬ 
cles) and kaons (empty blue squares). The measured mean 
values agree well with the calculated dependencies of eq. 
0, shown by the dashed and full lines, respectively. No 
momentum dependence has been found for the pion reso¬ 
lutions. For the kaons, we observed some fluctuations due 
to the smaller statistics of the data sample, however the 
fitted values are the same as for the pions within about 
Icr. 

The number of sigma separation between the pion and 
the kaon Cherenkov angle peaks can be computed using 
the standerd gaussian estimator 


TIq- 


Vc-Vc 

K + 0'^]/2‘ 


( 10 ) 


We obtain Ucr > 4 in the whole momentum range, cor¬ 
responding to a pion rejection factor of about 1:500, as 
required for the CLAS12 RICH. Non-gaussian tails in the 
Cherenkov angle distributions could reduce the real re¬ 
jection power. Part of the tails are due to non optimal 
tracking reconstruction (in particular events with more 
than one track). In addition, events were accepted start¬ 
ing from a minimum number of photoelectrons equal to 4. 
This very low number could be increased up to more real¬ 
istic values of 6 or 7 with a minor expense of less than 10% 
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Fig. 18. Top plot: Number of p.e. per Cherenkov ring mea¬ 
sured with the fourteen H8500C (red full line) and with the 
fourteen H8500C-03 (blue dashed line). Bottom plot: Cheren¬ 
kov angle resolution as a function of Npe measured with the 
fourteen H8500C (red squares) and with the fourteen H8500C- 
03 (blue triangles) MAPMTs. Data for 8 GeV/c momentum 
pions and aerogel with n — 1.05 and trad = 20 mm. 


of efficiency. The requirement of one track and at least 7 
p.e. provides a 1:500 rejection power as measured on the 
real distribution of tagged pions. 


5.3 UV photon resolution 


The array of MAPMTs includes 14 normal glass window 
H8500C and 14 UV-enhanced glass window H8500C-03, 
alternated along the ring. The effect of the UV photons 
can be estimated by comparing the photon yield and Che¬ 
renkov angle resolution using only the MAPMTs of the 
same type. In Fig. [181 we report this comparison for pi¬ 
ons with a momentum of 8 GeV/c. The top plot shows 
the number of p.e. per event. We see that the H8500C- 
03 MAPMTs detect about one hit per event more than 
the H8500C MAPMTs. The bottom plot shows the Che¬ 
renkov angle resolution as a function of Npe. In spite of 
seeing more light, the H8500C-03 provide a worse resolu¬ 
tion, because UV photons suffer more scattering inside the 
aerogel and in the air gap and also because of the larger 
chromatic dispersion (see Sect. lEH). Fitting the data with 
eq. (0, we extracted aipe = 4.01 ±0.01 mrad for the 
H8500C and aipe = 5.03 ± 0.02 mrad for the H8500C-03. 

Given these results, the H8500C MAPMTs have been 
chosen for the CLAS12 RICH application. 


5.4 Aerogel refractive index 


Measurements have been performed with aerogel tiles with 
nominal refractive index n = 1.04,1.05,1.06 and thickness 
trad = 20 mm. In Fig. na we show the Cherenkov angle 
measured for the three refractive indices (full circles) com¬ 
pared with the expected dependence (curve). Deviations 
from the expected value, represented by the curve, indi¬ 
cate a difference of the actual refractive index from the 
nominal one. For the two smaller refractive indices, the 
measured difference is below 0.1%, while for n = 1.06 a 
difference of the order of 0.2%was found, within (although 
close to) the typical tolerance achievable with the present 
production technology. 

In the top plot of Fig. [20] we show the measured Npe by 
the full circles, while the empty squares show the results 
rescaled by the effective ring coverage, which is higher for 
lower n. The number of p.e. emitted at a given wavelength 
and after a radiator thickness x can be computed by fold¬ 
ing the Cherenkov emission spectrum of eq. o with ex¬ 
perimental effects, as 


£n 

dxdX 


d?N 


dxdX 


QE{X)epee M>^Tcos(r,c) 


( 11 ) 


where QE{X) is the MAPMT quantum efficiency [19], 
Cpe = ^coiEpack^rmSPE IS the total p.c. detection efficiency 
and A{X) is the aerogel transmission length, obtained from 
the measured transparency T(A) using the relation 


T(A)=e“^. (12) 

Integrating equation m over A and x, we obtained the to¬ 
tal number of p.e. values shown as the blue dashed curve 
in the top plot of Fig. |20l which is in sufficiently good 
agreement with the measurements for the refractive in¬ 
dex n = 1.04,1.05. The measurement is significantly lower 
than the expectation for n = 1.06. This may indicate for 
this aerogel tile larger contributions from effects not con¬ 
sidered in eq. ([IT]) (for example surface quality or density 
variations). 

In the bottom plot of Fig. |20l the full red circles show 
the measured Cherenkov angle resolution cr^, while the 
empty blue squares show the single photon resolution aipe 
extracted from the Npe dependence using eq. (|8|). While 
the results for the lower refractive index tiles are compa¬ 
rable, worse resolution has been found for the aerogel with 
n = 1.06. The larger value of is the combination of the 
lower Npe and of the worse single photon resolution aipe. 

These results indicate the n = 1.05 as the optimal 
choice for the CLAS12 RICH, because it provides the 
highest tt/K separation power by combining the biggest 
number of p.e. with the best angular resolution. 


5.5 Chromatic dispersion 

The chromatic dispersion of the aerogel used in the tests 
has been studied by measuring the Cherenkov angle as a 
function of the photon wavelength. This has been done by 
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Fig. 19. Cherenkov angle measured with aerogel tiles of of 
thickness t = 20 mm and different nominal refractive index 
(full circles) compared with the expected values (dashed line). 
Data for pions with 8 GeV/c momentum. 




Fig. 20. Top plot: Average number of p.e. of the data (full cir¬ 
cles) and rescaled to full ring coverage (empty squares) and ex¬ 
pected Npe (dashed line). Bottom plot: Cherenkov angle resolu¬ 
tion (full red circles) and single photon resolution aipe (empty 
blue squares). Data for pions with 8 GeV/c momentum for 
aerogel radiators with same thickness trad — 20 mm and three 
different nominal refractive indices. 


placing, just after the aerogel radiator, optical filters [32], 
which select photons in a band of wavelengths. 

The Cherenkov angle reconstruction basically extends 
the eq. (HD to include refractions at the surfaces separating 
the air and the filters. At each wavelength, the refractive 
index is then obtained from the Cherenkov angle through 
the eq. ([7|). 

In Fig. mi we show the results for the aerogel with 
nominal refractive index n = 1.05 and trad = 20 mm 
obtained with 8 GeV/c momentum pions. The full circles 
are our data, with the vertical bars representing the la er- 



Fig. 21. Refractive index as a function of the Cherenkov light 
wavelength for the aerogel with n — 1.05 and trad = 20 mm 
obtained selecting pions with 8 GeV/c momentum. The full 
curve and the gray area represent the result of a fit with the 
dispersion relation eq. (O, whose fitted parameters are shown 
in the legend, and the la error band. The dashed pink curve 
is the result of the calculation described in the text. 


ror on the mean refractive index and the horizontal ones 
representing the wavelength interval of each measurement. 
The curve represents a fit using the Sellmeier parametriza- 
tion [33] of the wavelength dependence 

= l + (13) 

that reproduces the data well. The gray area in the plot 
represents the la error band of the fit. From this fit, we 
obtained a refractive index n/it (400 nm) = 1.0494T0.0004 
nm at the wavelength Aq = 400 nm. The variation in the 
wavelength range between 300 and 600 nm is about ±0.2% 
with respect to the value at 400 nm. 

According to the quotation provided by the Novosi¬ 
birsk factory, the refractive index at the wavelength Aq 
can be obtained by using the equation 

n^(400 nm) = 1 + 0.438p. (14) 

Using in the measured aerogel density p = 0.230 g/cm^, 
we got the value nquoted = 1.0492, in agreement with our 
measurements. 

Following Vorobiev and collaborators [3l], the A de¬ 
pendence of the aerogel refractive index can be computed 
from the well known values of its main components, the 
air and quartz, through the relation 

Tl{X) = AUairi^) A^Hguartzi^) - 

The normalization constant A can be computed, if the 
aerogel refractive index is known at the wavelength Aq , as 

_ '^quartzi^o) (16) 

"f^quartzi^o) '^air (Ao)' ^ ^ 

Using Ao = 400 nm and n(Ao) = n/it(400 nm), we ob¬ 
tained the result shown in Fig. [21] by the dashed pink 
line. It reproduced the measured wavelength dependence 
very well. 





















14 


S. Anefalos Pereira et aL: Test of the CLAS12 RICH prototype 


6 Analysis of the Cherenkov angle resolution 


The resolution of the measured emission angle of the Che¬ 
renkov photon can be written as 

^Ipe =<^l + ^d+<^l ( 17 ) 

with main contributions arising from the chromatic dis¬ 
persion (cTc) and from the uncertainties on the impact 
point measurement {(id) and on the emission point (ae). 
We neglect here other contributions like background pho¬ 
tons (largely removed in the analysis) or reconstruction 
issues (following the discussion on the Npe dependence of 
the measured SPE angular resolution). We also neglect 
the potentially important contributions associated to the 
aerogel optical properties (Rayleigh scattering, surface ef¬ 
fects, etc.) that are more difficult to determine and which 
require specific experimental studies. 

The chromatic term is due to the variation of the aero¬ 
gel refractive index with the photon wavelength and it can 
be further decomposed in a contribution (a*) due to the 
photon emission angle and a contribution (crj) due to the 
refraction at the crossing of the aerogel surface and can 
be written as 


= {<? + {<) 


r\2 



where the variation of the air refractive index has been 
neglected. The variation of the aerogel refractive index 
an is calculated by using the fit to our measurements, 
as shown in Fig. [211 weighted with the spectrum of the 
detected photons computed by integrating eq. dH) over 
the radiator thickness. 

The contribution due to the uncertainty on the impact 
point can be written as 


CTd = 



(19) 


where, assuming a uniform photon hit probability within 
the square MAPMT pixels, the Cherenkov hit radius reso¬ 
lution is aR = d = 6 mm being the MAPMT pixel 

size. 

The contribution due to the uncertainty on the emis¬ 
sion point X can be written as 


= 



( 20 ) 


where it was conservatively assumed that, on average, the 
photon has been emitted at half of the aerogel radiator 
thickness, thus ax ~ trad/V^- 

In Fig. 122 we show the result of the calculated aipe 
for the 14 H8500C with normal glass window (top plot) 
and for the 14 H8500C-03 with UV-enhanced glass win¬ 
dow (bottom plot). The various curves are plotted as a 
function of the gap length, for an aerogel radiator thick¬ 
ness trad = 20 mm. The emission point and the pixel size 
contributions have similar size, decreasing with the gap 



Fig. 22. Calculated contributions to the Cherenkov photon 
angular resolution aipe as a function of the air gap length for 
the fourteen H8500C with normal glass window (top plot) and 
for the fourteen H8500C-03 with UV-enhanced glass window 
(bottom plot): the curves represent the calculated resolutions, 
while the points represent our measurement. Data for 8 GeV/c 
pions and aerogel with n = 1.05 and trad = 20 mm. 


length, and give a significant contribution only for a gap 
length of the order of 1 m or less. They are also indepen¬ 
dent of the MAPMT type. The chromatic term is almost 
fiat and is the dominant one, except for small gap length 
below approximately 1 m. 

We also show in the figure the measured total aipe 
values at the gap length tair = 994 mm of our prototype. 
The calculation results for the H8500C MAPMTs agree 
with the measurements within less than 10%. The total 
of possible missing contributions can be at most of the 
same order of magnitude of the smaller ones in the top 
plot of Fig.|22l For the H8500C-03 data, the discrepancy is 
slightly larger, even though below 15%, indicating that the 
missing contributions (as for example the Rayleigh term) 
might be larger in this case. The almost 30% increase in 
the chromatic term due to the larger wavelength interval 
extending into the UV region is anyway not compensated 
by the about 20% increase in the photon flux measured 
for the H8500C-03 MAPMTs. This supports the results of 
Sect. 15.31 that the UV photons may worsen the resolution 
in spite of increasing the number of photoelectrons. 


7 Conclusions 

A RICH detector has been proposed to extend the PID 
capabilities of the CLAS12 spectrometer of the Hall B at 
JLab in the momentum range between 3 and 8 GeV/c. 
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The detector presents an innovative hybrid configuration 
with an aerogel radiator, MAPMTs as photon detectors 
and mirrors, in which the Cherenkov rings are imaged ei¬ 
ther directly (for forward particles) or after two mirror re¬ 
flections and a double pass through the aerogel (for large 
angle particles). 

We reported here on the beam tests of a large scale 
prototype of the detector performed at the T9 hadron 
beam line at CERN for the proximity focusing configu¬ 
ration. The setup used multianode photomultiplier tubes 
Hamamatsu H8500 (a device not optimized for the sin¬ 
gle photon regime) and readout electronics based on the 
MAROC3 chip. The results showed that they provide good 
performances in the single photon regime. The radiator 
was aerogel with various refractive index values. 

In the most challenging conditions expected for the 
CLAS12 RICH, i.e. 8 GeV/c momentum pions, we mea¬ 
sured Npe = 12.02 ± 0.02 p.e. per ring and aipe = 4.58 ± 
0.02 mrad. These results, that provided the required pion 
to kaon rejection factor of 1:500, have been obtained with 
an aerogel radiator with n = 1.05 and trad = 20 mm, the 
final choice for the CLAS12 RICH. The measured resolu¬ 
tions are in reasonable agreement with the expectations. 

The tests also included measurements of several as¬ 
pects of the foreseen RICH detector: contribution of UV 
photons to the angular resolution, aerogel refractive index 
dependence and chromatic dispersion. These measurement 
showed that all these effects are effectively under con¬ 
trol, without affecting the required performances of the 
CLAS12 RICH. 
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